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Abstract 

This study presents a semi-analytical solution 
method to analyze the geometrically nonlinear response 
of bonded composite single-lap joints with tapered 
adherend edges under uniaxial tension. The solution 
method provides the transverse shear and normal 
stresses in the adhesive and in-plane stress resultants 
and bending moments in the adherends. The method 
utilizes the principle of virtual work in conjunction with 
von Karman’s nonlinear plate theory to model the 
adherends and the shear lag model to represent the 
kinematics of the thin adhesive layer between the 
adherends. Furthermore, the method accounts for the 
bilinear elastic material behavior of the adhesive while 
maintaining a linear stress-strain relationship in the 
adherends. In order to account for the stiffness changes 
due to thickness variation of the adherends along the 
tapered edges, their in-plane and bending stiffness 
matrices are varied as a function of thickness along the 
tapered region. The combination of these complexities 
results in a system of nonlinear governing equilibrium 
equations. This approach represents a computationally 
efficient alternative to finite element method. 
Comparisons are made with corresponding results 
obtained from finite-element analysis. The results 
confirm the validity of the solution method. The 
numerical results present the effects of taper angle, 
adherend overlap length, and the bilinear adhesive 
material on the stress fields in the adherends, as well as 
the adhesive, of a single-lap joint. 
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Nomenclature 


A (p) 

= area of adherends and adhesive 

p (p) 

= boundary of adherends and adhesive 

p(p) 

m 

- m th boundary segment 

n (P) 

11 m 

= unit normal to the m th boundary segment 

U 

= reference planes 

h (p) 

= half thickness of adhesive and adherends 

£(P) 

= coordinates on reference planes 

ti p) 

= ply thickness 

e[ p) 

= ply orientation 

E[ p) , E ( t p) , G[ p) , = ply material properties 

o (p) 

= coefficients of reduced ply stiffness 


matrix 

A <p) , D (p) 

, B ip) = in-plane, bending, and coupled 


stiffness matrices 

f( a ) 

eff 

= effective adhesive Y oung ’ s modulus 

Q(a) 

= effective adhesive shear modulus 

v (a) 

= adhesive Poisson’s ratio 

T lff > r$ 

= effective transverse shear stress and shear 


strains 

K 

= degree of the B-spline functions 

f{p) 

V 

= pre-selected knot points of B-spline 


functions 

T m (x; t ip) , K) , T n ( y ; t (p) , K) =B-spline functions 

Mf 

r (P) 

U a(mn) 

= extent of of B-spline functions 
= unknown coefficients in B-spline 

functions 

C (P) 

= unknown vector of coefficients, c ( a p J n 

C ? ty ? C 

= external in-plane tractions 

m x , m y 

- external bending moments 

N (/?) M (/?) 

= in-plane stress resultants and bending 


moments 

u ( /\ u ( y p) 

, U\ p) =global displacements 

u (p) 

m 

= prescribed displacements and slopes 

q ip) 

= unknowns in governing equations 

Jpl 

e afi 

= strain components 

e (p) 

= total strain vector 

Up) Up) 

K xx 9 K yy 

, k[ p) = bending strain (curvature) resultants 
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£ (p) 

, £ ip) 

XX 

’ yy 

~(p) 

K (p) 


“ ZZ 

9 °ZZ 

cr (a) 


y (a) 

/ az •> 

“zz 

c ( a ) 



A (p) A' (p) = 

1 am 9 1 zm 

A (p) 

m 

sw e = 


sw™ 


8W c (p) 

K 


f = in-plane strain resultants 

in-plane strain resultant and curvature 

vectors 

adhesive normal stress and strain 
adhesive transverse shear stresses 
adhesive transverse shear and normal 
strains 

transverse shear and normal strains in 
adhesive 

Lagrange multiplier functions 
matrix of Lagrange multiplier functions 
virtual work of internal and external 
forces 

internal virtual work of adherends and 
adhesive 

virtual work of boundary reaction forces 
potential energy of the constraint forces 


Introduction 

The reduction of the transverse shear and normal 
stress concentrations along the edges of adhesive 
bondlines is important in order to prevent premature 
failure of the bonded joint. The determination of the 
complete stress and strain fields in composite bonded 
single-lap joints presents difficulties arising from the 
step-wise geometry, material property variations, 
laminated construction of the adherends, the bilinear 
material behavior of the adhesive, as well as the effect 
of tension induced stiffening (geometrically nonlinear 
effect) on the bending deformation of the adherends 
subjected to uniaxial tension. Also, the local stress 
variations near the ends of the overlap region are 
characterized by very high gradients or even 
analytically predicted singularities. The sharp gradients 
of the stress components depend on the elastic 
properties of the adherends and adhesive as well as the 
joint geometry. These peak transverse normal and 
shear stresses in the adhesive can be reduced by 
tapering the adherends toward the ends. 

In order to facilitate the use of lap joints in present 
and future structures, the analysis of the geometrically 
nonlinear response of bonded single-lap joints has 
received considerable attention over the past two 
decades. Previous analyses of bonded lap joints can be 
categorized as “shear-lag” and “finite-element” models. 
An extensive review and in-depth discussion of the pre- 
vious investigations can be found in articles by Tsai and 
Morton, 1 Ding and Kumosa, 2 and Osnes and 
Andersen. 3 Due to the aforementioned complexities, the 
majority of the investigations have utilized the finite 
element method in determining the stress and the strain 
field in a bonded lap joint. However, in many of these 


investigations, the three-dimensional description of the 
bonded lap joint was simplified to a two-dimensional 
analysis under certain assumptions. 4 " 7 

In order to enhance computational efficiency, 
Penado 8 introduced an approach based on the sub- 
structuring technique. In this approach, the general 
response of the bonded lap joint is obtained analytically 
from the solution of force-moment equilibrium 
conditions. The analytically evaluated force-moment 
values at the overlap ends are then used as the natural 
boundary condition for a highly detailed two- 
dimensional finite element analysis of the overlap joint 
under the assumption that the overlap ends are simply 
supported. 

Combining the shear-lag model of Goland and 
Reissner 9 with a detailed finite element modeling of the 
adherends with three-dimensional elements, Edlund and 
Klarbring 10 employed the principle of virtual work to 
analyze the geometrically nonlinear response of bonded 
single-lap joints with a linearly elastic adhesive. The 
shear-lag model approximates the transverse shear and 
normal strain components in terms of the relative 
displacements of the adherends. An alternative to the 
shear-lag model is to model the adhesive with one or 
two layers of brick elements and assemble these 
elements with the brick elements of the adherends. 
However, as mentioned before, the size of the brick 
elements used in the adhesive might introduce aspect 
ratio problems. The recent investigations by Pandey 
and Narasimhan 11 and Narasimhan and Pandey 12 
utilized this approach for solving the three-dimensional 
large deflection analysis of single-lap joints with 
viscoelastic adhesive behavior. 

In the finite element analysis, the adhesive requires a 
highly refined mesh in order to keep the proper aspect 
ratio between the elements in the adherends and 
adhesive. Therefore, the major advantage of the two- 
dimensional finite element models over the three- 
dimensional models is the significant reduction of the 
number of degrees of freedom. 

Furthermore, for an incremental-iterative solution of the 
governing equations in which the global stiffness 
matrix is repeatedly calculated, the three-dimensional 
finite element analysis of the entire domain becomes 
computationally demanding. 

Thus, it is beneficial to have an efficient special- 
purpose analysis method that can be used to conduct 
extensive parametric studies in a timely manner and at 
relatively low computational costs. However, there is 
no analytical or semi-analytical approach for 
determining the three-dimensional response of the 
geometrically nonlinear analysis of bonded single-lap 
joints. Therefore, the goal of the present study is to 
develop a three-dimensional analysis method that is 
well suited for parametric studies that accurately 
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predicts the geometrically nonlinear behavior of a 
bonded single-lap, tapered, composite joint subjected to 
uniaxial tension. In particular, this study focuses on the 
effects of geometric nonlinearity, tapering of the 
adherend thickness, changes in the overlap length of the 
adherends, and the linear and bilinear elastic behaviors 
of the adhesive on the in-plane stresses in the adherends 
and the transverse normal and transverse shear stresses 
in the adhesive of the joint. 

In the remainder of this paper, details of the analysis 
method are presented and results from this approach are 
discussed. First, the boundary- value problem is 
defined. Next, the analysis details and numerical 
solution procedure are described. Then, results for two 
bonded lap-joint configurations are presented. 

Problem Statement 

The bonded single-lap joint configuration consists of 
two rectangular composite adherends bonded through a 
thin layer of adhesive, as shown in Fig. 1. While the 
adhesive has a uniform thickness, the composite 
adherends have tapered edges. The tapered edges of the 
adherends are used to reduce stress concentrations 
along theedges. Both the geometry and kinematics of 
the adherends and adhesive are described by utilizing 
the global coordinate system (x, y, z), as shown in Fig. 
1 . 



Fig. 1 Description of geometry, fiber angle, and 
loading of a bonded single-lap joint subjected to 
uniaxial loading. 

The lower and upper adherends are identified by the 
super scripts u and i , and their associated area and 
boundaries are denoted by A (u) and A (£) and T (M) and 
T (0 , respectively. The adhesive is identified by the 
subscript a and its associated area and boundaries are 
represented by A (a) and T (a) , respectively. As shown 
in Fig. 1, both the adherends and adhesive have a 
rectangular boundary geometry. Hence, their 
boundaries can be represented by dividing their entire 
boundaries into four straight boundary segments, i.e, 


r (p) =Z r lf ) (p = u,l,a) ( 1 ) 

m = 1 

in which T {p) denotes the m th boundary segment 
around the adherend and adhesive boundary. Along the 
adherend and adhesive boundaries, T p (p = u,£,a ) , 
the unit normal to the m th boundary segment is 
represented by n^ } , with components n { x p J and n$ in 
the x- and y-directions, respectively. The unit normal, 
n^f } , makes an angle, </> ( m p) , with respect to the positive 
x-axis, as shown in Fig. 1. 

The exterior edges of the adherends are subjected to 
both in-plane tractions and bending moments. The in- 
plane external tractions include components t x , t y and 

t z , and the external bending tractions include compo- 
nents m x and m y . The traction components are defined 

with respect to the (x, y, z) structural coordinates, and 
their positive-valued directions are shown in Fig. 1. In 
the adherends and adhesive, the global displacement 
components in the x-, y-, and z-directions are denoted 
by U x p) , U y p) , and U ( z p) , with p = u 9 £, a , respectively. 

The laminated adherends are made of specially 
orthotropic layers. Each layer has a thickness of t[ p) and 

orientation angle of 0 { k p) {p-u,£), which is defined 
with respect to the positive x-axis, as shown in Fig. 1. 
Also, the orthotropic material properties of each layer 
include the elastic moduli E[ p) and E P p) , shear modu- 
lus G { l p) , and Poisson’s ratio , where L and T are 
the longitudinal and transverse directions, respectively. 



Fig. 2 Bilinear behavior of the adhesive in terms of 
effective transverse shear stress and transverse 
shear strain. 
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The adhesive material is isotropic, homogeneous, and 
elastic with a bilinear relation between the effective 
transverse shear stress, , and effective transverse 

shear strain, y { ^ , as shown in Fig. 2. The effective 
transverse shear stress and strain are defined by 


=V Cr xz )2 +cr yf 2 

(2a) 

riff =p ( x f + r ( y f 

(2b) 


in which and <j ( y a J represent the components of the 

transverse shear stress and y^J and y^ represent the 

components of the transverse shear strain in the 
adhesive. As shown in Fig. 2, the initial shear modulus 

of the bilinear adhesive behavior is denoted by , 
and it reduces to G[ a) after the characteristic transverse 
shear strain, y { c a) . It has a Poisson’s ratio of v (a) . With 
these parameters, the bilinear relationship between the 
effective transverse shear stress, , and effective 

transverse shear strain, y ^ , can be expressed as 



where H(y { e ^ -y ( c a) ) is the Heaviside step function. 

In accordance with this relationship, the transverse 
shear stresses, cr^J , and strains, y^J , are related by 


Eg=2Gg(l + v M ) (7) 

While the thickness of the adhesive is uniform and 
denoted by 2 h (a) , the adherends have variable 
thicknesses, h ip) (x,y ), with p = u,£, due to their 

tapered shape near the edges. The thicknesses of the 
adherends are defined as discrete linear functions of the 
in-plane coordinates. The reference planes of the 
adherends and adhesive, denoted by z ip) , 

with p = u,£,a , are shown in Fig. 3. 

The problem posed here concerns the development 
of a three-dimensional semi-analytical method to 
determine the displacement and stress fields in bonded 
single-lap joints while including the effects of 
geometric nonlinearity and bilinear elastic adhesive 
material behavior. The capability of the method is 
demonstrated by considering first a quasi-isotropic 
tapered adherends and a bilinear adhesive material 
behavior. The second configuration concerns the effect 
of the overlap length of the adherends on the general 
response of both the adherends and adhesive. The lap 
joints are simply supported along the left and right 
edges of the upper and lower adherends, respectively. 
The left edge of the upper adherend is also restrained 
against any horizontal movement. Furthermore, the 
mid points of the left and right edges of the lower and 
upper adherends are restrained to move in the y- 
direction so that the rigid-body movement of the lap 
joint is completely suppressed. 

Solution Method 

The present three-dimensional geometrically 
nonlinear analysis method is based on the principle of 
virtual work. The displacement components are 
approximated in terms of the B-spline functions 13 in a 
double series representation as 


( a = x,y ) (4) 

in which the parameters G ^ represent the effective 
shear modulus of the adhesive, defined as 


m ( p) My 

< P) = Z Z #L) T m T n (y,ty,K) (8) 

m = 0 n = 0 

in which c ( a p) mn) , with p = u,£ and a = x,y,z , are the 


G (a) = 



( 5 ) 


Furthermore, the transverse normal stress, cj { z f , and 
strain, , in the adhesive are related by 


r («) 


. (a) 

' e// b zz 


( 6 ) 


in which is the effective Young’s modulus 

expressed as 

4 


unknown coefficients. The parameter M ip \ with 
P = x,y and p = u,l , specifies the extent of the series 
and the knot vector t {p) contains pre-selected knot 
points (coordinates) in the direction of P . These knot 
points are used to increase the accuracy of the B-splines 
at certain locations. Also, the parameter K controls the 
degree of the B-spline functions and, consequently, the 
continuity of the field variable. The K th -order B-spline 
function is comprised of the (K - l) th -order of 
polynomials. 
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Fig. 3 The reference surface and kinematics of the bonded single-lap joint. 


The details of the B-spline functions along with the 
definitions of the parameters t { p and K are given in 

the Appendix .These displacement functions can be 
expressed in matrix form as 

u[ p) 

u (p) =y (P)\(P) (9 a, b, c) 

U (P) =Y (p)\(p) 

in which the vectors c^ , with a = x,y,z , contain the 
unknown coefficients (generalized coordinates) c ( a p J n . 
The known vectors, V (/7) , are expressed in terms of the 
B-spline functions T m (x;t[ p) ,K) and T n (y;t (p \ K) . In 

matrix form, the approximate displacement 
representations of Eqs. (9a-c) are rewritten as 

u ip) =\ { x p)T q ip) 

u (p) =y(p) T q (p) ( 10 ) 

u[ p) =\ { z p)T q ip) 

in which the known vectors, \ { a p) , with a = x, y, z , are 


defined as 

vf )T ={v (/0 ,o,o} 

(11a) 

v‘ p)T ={o,v <p) ,o} 

(lib) 

\y T =|o,o,v <,,) | 

(11c) 

The unknown vector q ip) is defined as 

q,,.’ = { C WV«V«’} 

(12) 


Note that the series representation of the displacement 
components is not required to satisfy any type of 
kinematic admissibility. 

Displacement Components 

The adherends interacting through the adhesive, 
which sustains transverse normal and shear 
deformations but not in-plane deformation, are 
subjected to in-plane and bending deformations but not 
transverse normal and shear deformations. The 
transverse normal and shear strain components in the 
adherends are disregarded because they are thin. 
Therefore, the in-plane strain components in the 
adhesive, and the transverse normal and shear strain 
components in the adherends, are not included in the 
derivation of the kinematic relations. 

In accordance with the Kirchhoff plate theory, the 
global displacement components, U {p) , U {p) , and U[ p) 
in each of the adherends are defined as 


uy ( x,y,z ) = u ( y (x,y)-C P) h M u% 

(13a) 

U ( z p) ( x,y,z) = u[ p) (x,y ) 

(13b) 


for which p-u,t and a - x,y , and the displacement 
components, U [ p) , U ( y p) , and u[ p) , are defined on the 
reference surfaces with respect to the global Cartesian 
coordinates ( x ip \ y {p \ z ip) ) , as shown in Fig. 3. In Eq. 

(14), the subscript after a comma indicates 
differentiation with respect to the variable. The 
coordinate p p) located on each of the reference planes 
is defined as 

z-z {p) 

S (P) =— h ur (p = u,£) (14) 


5 

American Institute of Aeronautics and Astronautics 


and varies in the range -1 < £ ip) < 1 , with p -u, £ . 
The thicknesses of the adherends are specified by 
2 h (p \x,y) , in which h ip \x,y) = h^ p) in the untapered 
sections of the adherends. Also, the location of the ref- 
erence planes with respect to the global coordinate sys- 
tem (x,y,z) are defined by z (p) , with p = u,l , which 
are located at the mid-surfaces (with respect to the unta- 
pered thickness) of the adherends, as shown in Fig. 3. 

Strain-Displacement Relations 

The strain measure for the adherends is based on the 
modified form of Green’s nonlinear strain displacement 
relations in conjunction with von Karman assumptions 
for large deformation of plates 14 . Therefore, the strain 
components in the adherends, e {p ^ (p = u,l; a,P~ 
x, y ), can be expressed as 


( P ) = £ { P ) _ r^h {p) K (p) +-U P) f 

xx XX *5 xx ^ \ z ’ x ) 

(15a) 

e (p) =s ip) -C (p) h {p) tc (p) +-(u (p} ) 2 

yy yy * yy 2 ^ z,y ' 

(15b) 

e ( p ) = L y (p) +i £(p'>h ip) K (p) +-u (p) u (p) 

xy ^ ' x y 2~^ xy 2 z,x z,p 

(15c) 


nonlinear part of the in-plane strain resultant vector, 
£ {p) , in the form 



Although the bending deformations (curvatures) are 
only linearly related to the out-of-plane displacement 
component, u ( z p) , for consistency, a zero-valued vector 
is employed to represent the nonlinear part of the 
curvature vector, , as 

k ( ; } ={o,o,o} (i9) 

Substituting for the derivatives of the displacement 
components from Eq. (10), the linear and nonlinear 
parts of the in-plane strain resultant and curvature 
vectors can be expressed as 

e[ p) = Ii£q ip) , k[ p) = (20a, b) 

e^=L^«X, <=0 (20c, d) 

where 


in which 


C-(P) 

^XX 

= u (p) 

x,x 

trip) - ip) 

’ XX U z,xx 

(16a, b) 

~(P) 

b yy 

= u (p) 

y.y 

^.(P) _ _„(P) 

9 K yy ~ u ^yy 

(16c, d) 

y(P) 
/ xy 

= u ip) 

x,y 

+ u ip) ■ K (p) - -7u (p) 

+ U y,X ’ K xy - ZU z,xy 

(16e, f) 


where s [ p) , £ y p) , and y { x p) represent the in-plane strain 
resultants and /c x p) , K y p) , and k[ p) represent the 

bending strain (curvature) resultants on the reference 
surfaces. Also, the in-plane and bending (curvature) 
strain resultants constitute the components of the linear 
part of the in-plane strain, z {p) , and curvature, , 
vectors in the form 


y(/T 

x,x 


L (p) = 


y(p) T 

y,y 

y(^) T + y (^) T 

J x,y 1 y,x 


y(/T 
T z,xx 


r (p) — 


y(^) T 

z>yy 

2V (/7)T 

^ Y Z,Xy 



( p ) y (p) t 

z,x z,x 

U (P) y (pf 
z,y z,y 

u (p)y( p ) T +u (p)y( p ) T 

z,x z,y z,y z,x 


(21a) 


(21b) 


(21c) 


The vectors of strain resultants defined in Eq. (20) can 
be combined in a compact form as 







( p ) 

y,y 


,u 


(p) 

x,y 


(p) ,-u (p) 

z,xx ? z,yy 


+ «#} 
-2<} 


(17a) 

(17b) 


Similarly, the nonlinear terms appearing in the strain 
components, e iP p (a, P = x, y ) , are included in the 


z {p) =li p) q ip) ( a = L,N ) 
where 


e 


t p) T 



ApY 


)• 



( 22 ) 


(23a, b) 
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(23c) 


in which 


L ( ^(q (/?) ) = 


0 


Furthermore, the linear and nonlinear parts of the strain 
vectors, and e^ } , can be added to form the total 
strain vector as 



v^ T 

vf )T 

o r 


(p 


£,u) 


(30a) 


e (p) = e[ p) + e (p) = [ii/ 0 + ll p) (q (/,) )] q <p) 
= H <P ) (qWjqd’) 

where 


(24) 



h (p) \ { z f 

h (p) \ ( z f 

(_l)Vy(^ )r 


(p = £,u) 


(30b) 


H (/,) (u [ p) ) = li p) + (q (;,) ) (25) 

For the adhesive bonding of the two adherends, the 
displacement components are assumed to vary linearly 
through the thickness 1 . Although the adhesive 
undergoes the same magnitudes of the in-plane and 
transverse displacements as those of the adherends, the 
strain measure is based on a linear shear-lag model, 
where the transverse shear strain and the normal strain 
components in the adhesive are expressed as 


=^[U?(x,y,0)-U?(x,y,0)] 


+ 


+ 


2l 


|l + C ( " l )te',;|r,v,0] a = (26a) 


(l-r’)C'2(a,J',0( 


(26b) 


Finally, substituting from Eq. (10) for the displacement 
components in Eq. (26) leads to the strain vector 
containing the transverse shear and normal strain 
components in the adhesive as 

£ («) =L ( : ) q (,)_ L (0 q (0 (27) 
where 


Stress-Strain Relations 

The external uniaxial in-plane loads acting along the 
boundary of the adherends, as shown in Fig. 1, result in 
not only in-plane stresses but also in bending moments 
in the adherends due to the couple formed by the 
eccentric location of the applied in-plane forces. The 
resulting transverse normal (peeling) stresses created in 
the adhesive are also due to this couple, which causes 
the bending deformations in the adherends. 

The in-plane stress resultants and bending moments 
generated by the applied external loading in the 
adherends are related to the in-plane strain resultants 
and curvatures, which are defined on the mid-surfaces 
of the upper and lower adherends through the 
constitutive relation as 

f N (?) i 

|m ( ^J 

where 


X {p) (x, y ) B (p) (x, y) 
B (p \x,y) D (/,) (x,j;) 



iv p, i 

r 

l K<P) J 


(31) 


4 P) ( x >y) = h(p) -ti p) ) Q% 

k=\ 

4 P) ( x >y) = \{ h(p) (M) 2 Ife 2 - 4 P)1 ) (32) 

D!j p) (x,y) = l -(h^ (x,y)j f^ 3 -ft * ) 

with 


-{/is?,*?} 

and the matrices li p) are defined as 


(28) 


(29a) 

(29b) 


MV = z _4^y)z41 

k h (t) (x,y) 

(k = l,...,N e ;z (t) -h (t) <z k <z (t) + h (<) ) (33a) 

mu) = fA^y)z41 

^ h M (x,y) 

(k = 1, . . . , N u ; z (u) - h (u) < z k < z (u) + h (u) ) 


T (m) — T 0) , T (w) 
^ as 

tW = tW_tW 

ax 
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(33b) 



(38b) 


In Eq. (31), the matrices A (/?) , D (/?) , and B (p \ with 
p = u,£, are associated with in-plane, bending, and 
coupled in-plane and bending behaviors of the 
adherends, and Q\^ k) (p-u,£) are the coefficients of 

the reduced stiffness matrix of the k th ply defined in 
the global (x-y) coordinate system. Note that the tapered 
adherend thickness, h ip) (p = u,l), varies as a function 
of the (x-y) coordinates. Hence, the material property 
matrices associated with the adherends, A (/?) , D (/?) , and 
B (/?) , are dependent on the in-plane coordinates. 

Furthermore, the ratio of the ply thickness to the 
adherend thickness is assumed to be constant, i.e., 
t ( k p \x,y)/ h {p) (x,y) = T k ip) = constant . In this case, the 

material property matrices, A ip) , B (p) , and B (/?) , 
become dependent only on the adherend thickness, 
h ip) (x,y) , with p = u,£. 

The relation given in Eq. (31) can be compacted in 
the form 




,T ={r { : ) y; ) y?} 

E ( " , (q w ,q < '’ , ) = 


C(q w >q w ) o 

0 G$( q w ,q w ) 

0 0 


E eff W ’<1 j. 


(38c) 


in which the expressions for (q (5) , q ip) ) = G ^ (y ( e ^ ) 

and E { ^ (q (,v) , q (p) ) = (y ^ ) are defined in Eqs. (5) 

and (7), respectively. 

Substituting for the expression for c (a) from Eq. (27) 
permits the stress-strain relations given in Eqs. (38) to 
be expressed in terms of the unknowns of the adherends 
as 


s (/>) =E (/>) e 0>) (p = u,l) (34) 

in which s (/7) , E (/?) , and e ip) are defined as 

s (rt T = j N (rt T 5 M (rtT } (35a) 

, A (P) 

E w = M “ (35b) 

g(P) J)0) V 7 


e (p)T = U (p) \ k ( p)T 


(35c) 

With the representation of e (p) in Eq. (24), the stress- 
strain relations given in Eq. (35) can be rewritten as 


s ( /h = E (p) H (p) (q (p) )q (p) (p = uj) 


(36) 


Because the adhesive does not sustain in-plane 
deformation, the in-plane stress components, <j\ 


O) 


<j ^ , and cf , are disregarded. The transverse shear 


A a ) and rr (a ) 

xy J 

stresses, cr^f and cr^f , and the transverse normal 

stress, (7 ( z f , are related to the corresponding strain 
components through a bilinear relation as 


(a) _ j£(a) c (a) 


s v 7 = 


where 


«,(«) 


= h (a) { 


(7 (a) ,(7 {a) ,(7 {a> \ 

xz ’ yz ’ zz ) 


(37) 


(38a) 


|(‘> = E (a) (q < “ ) ,q (,) )(L ( a “V“ ) -L ( «q (/) ) 


(39) 


Boundary Conditions 

Along the m th segment of the boundary of the 
adherends, denoted by T ip) , with p-u^, as shown in 
Fig. 1, the prescribed displacement components normal 
and tangent to the boundary, 5 9 and im) u { z p) ? 

and the slope normal to the boundary, {m) u[ P l , can be 
imposed as 


U (P) = ("OjjOO 
n n 

(p) _ (m) -(p) 


onr L (p = u,t;m = W,4) (40) 


U (P) = 

U (P) = (*)£(/» 
z,n z.n 


Utilizing the vector representations of the displace- 
ment components given by Eq. (10), these prescribed 
displacements can be expressed in vector form as 

y ( m p)T q (p) -U (p) =0 (p = u,£; m = 1,2, 3, 4) (41) 

where the matrix V 1 / 1 and the vector iff' are defined 
as 
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cos 

+ sin $„ v f’ )T 


-sin^ m 


y(p) T 

T X 


y(^) T = 


+ cos^V^ )T 


0 


0 


0 

yW T 


(42) 


0 


cos i,V^ )T 

+ sin^V«f 


where 7.(^;t ( r p) ,^) represents the j th term of the 
(if-l) th -degree B-spline functions defined along the 
straight boundary segment Y {p) , with t ( r p) being the 
associated knot vector. Also, ^ j(am) , with a = n,t,z , 
and /L y ' (zm) are the unknown Lagrange multipliers 
associated with each B-spline function, 7.(^;t ( r p) ,^r) . 

Substituting the expressions for the Lagrange 
multiplier functions from Eq. (46) into Eq. (45) and 
rearranging the terms, the constraint equations 
representing the prescribed displacements can be 
rewritten as 


;Cp) t 








{m) u {p) \ 

z,n ) 


(43) 


The boundary conditions in Eq. (41) are enforced as 
constraint conditions by introducing Lagrange 
multiplier functions, (£) , with a = n,t,z , and 

A ' z ip) ( t ) , defined along the m th boundary segment. 
These boundary conditions are written in integral form 
as 


\ A^(t){\^ T q^-u^}dt = 0 

p( P) 

L m 


(p = u,£; m - 1,2, 3, 4) 


(44) 


where the matrix A ( ^ } contains the Lagrange multiplier 
functions in the form 


A^(0 


a ( ^(o= 


Mm (0 


Mm (0 


A ’L p) (t) 


( p = m = 1,2, 3, 4) 


(45) 


The (unknown) Lagrange multiplier functions , A^ (/) , 
with a = n,t,z, and are assumed in 

polynomial forms as 


M p)T (c ( m p W p) -f^)=o (47) 

where 

M P)T ={M P M M P J\ 4f} (48) 

with 

={Mty Mty Mty «)} ( 49 ) 



with 


f, ( £> - f Tj(^ m ,K)tifdT (53) 

pCp) 

for (p = u,£;m = 1,2, 3, 4) . 

The constraint equations in Eq. (48) can be assembled 
to form a single matrix equation combining all of the 
constraint equations as 


(A2(f«>),Ar(<«») 

7=0 

Q? = w,4m = l,2,3,4) (46) 


^(? )T (c (/7) q (/,) - 

where 



f c (/' ) ) = ° (p = u,t) 

(54) 

M p) \ ■ 

4' ,T ) 

(55a) 

c</> T • 

.. Ct'> T l 

(55b) 
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4”’ =jfJf ,T > 4 P ’ T . -> 4 )T } (55c) 

The system of constraint equations in Eq. (54) is 
unique, provided the rank of the matrix C ip) is equal to 
the total number of constraint equations. Also, Eq. (54) 
can be treated as the potential energy of the reaction 
forces producing zero energy since C (/?) q (p) = 0 , and it 
can be referred to as the potential energy of the con- 
straint forces, V c , in the form 

K = Z ^ ( C (p) q (p) -f c (p) ) = 0 (56) 

p=u,t 

Governing Equations 

The governing equations are derived based on the 
principle of virtual work 

8W t = 8W e (57) 


where 

*H. (p \q (p) ) = l} p) +2li*\q ip) ) (62b) 

with p = u,£ . The external virtual work is expressed as 
the sum of the virtual work due to externally applied 
forces, 8W e (p) , and that arising from the boundary 

reaction forces, SW c (p) , i.e., 

SW e = £ 8W e (p) + SW c (p) (63) 

p=u,i 

The virtual work due to externally applied forces, 
SW e ip) , can be expressed in matrix notation as 

SW e (p) = Sq (p)T p ip) (64) 

in which 


where 8W i and 8W e represent the virtual work due to 

internal and external forces, respectively, of the bonded 
single-lap joint. 

The internal virtual work, SW i , is the sum of the 

internal virtual work of the adherends and the adhesive, 
i.e., 

§w t = SW i (u) + 8W™ + 5W t <a) (58) 

where the internal virtual work in the adherends and 
adhesive are expressed as 

SW; p) = J Se <p> \ ,p> dA = J Se (p)J E (p) e (p) dA (59) 

A p A p 

SW; a> = J Sz (a)T s (a) dA = J ^£ <a)T E (a) E <a) ^ (60) 

A a A a 


where A denotes the areas of the adherends 

(p = u,£) and adhesive (p = a). 

Substituting from Eqs. (24) and (27) and with the 
property of S[L ip) (q (/?) )]q (/?) = L ( ^ ) (q (/?) )^q (/?) , the total 

virtual strain vectors, 8z {a) and , are obtained as 


8* (a) =8q (u)T ll?- 

and 


Se (p) = [L ( £ p) + 2 L (p) (q ^ p) )] 8q (p) 
= *H <p) (q <p) )(?q </ ’ ) 


(61) 


(62a) 


P (p)T ={p ( / )T ,Pl P)T } (65) 

with 

pi p) = (66a) 

p(P) = J { yM Pz + vgm x + V%m y } dT (66b) 

p(p) 

The virtual work due to the boundary reaction forces, 
SW c (p) , is identical to the first variation of the potential 
energy expression in Eq. (56) as 


SW C (P) = -8V c 

= -SA (p)T (c (p) q (p) - f ( c p) ) (67) 

-Sq (p)T C (P)T A, ip) 


The virtual work due to the boundary reactions 
(constraint conditions) can be interpreted as the virtual 
work of the constraint forces, , over the virtual 

displacements, £(C (/7) q (/?) ) , of the adherends and the 
virtual work of the constrained displacements 
(boundary conditions), C (;?) q (;?) ■> over virtual 

constraint forces, 8A ^ P) . Although the term inside the 

parentheses in Eq. (67) is identical to zero, it is included 
in the virtual work expression in order to obtain a 
complete set of equations that contains both equilibrium 
equations and constraint conditions (kinematic 
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boundary conditions) along the boundary of the 
adherends. 

Substituting from Eqs. (59), (60), (64) and (67) 
while invoking the strain vectors, (p = u,l) and 
, from Eqs. (24) and (27) , and their virtual forms, 
8 ^ (p = u,l) and Sz^ a \ from Eqs. (61) and (62), 
and rearranging the terms, the virtual work expression 
of Eq. (57) can be rewritten as 

Sq (e)T (K w (q w ) + K;; ) (q w ,q (M) ))q (0 
+ <?q ( “ ,T (K <u) (q w ) + Kj> (q (f) , q , " ) )) q'" 1 
-<yq ( " )T K«(q ,,) ,q , " , )q , ' ) 

- Sq (l)T K (q <0 , q < “ ) )q < " ) (68) 

= Z Sq (p)T p (p) -SA (P> ' C (p) q (p) 

P =u,e 

-Sq (p)T C (p)T A (p) +SA (p)T C } 


where 

K (p) (q (p) ) = J *H wT (q <rt )E <p) H w (q w )^ 


(p = u,l) 


(69a) 


K$(q (, °,q (0 ) = J L^E^q^q^JL^ 

A a 

= 1,2) (69b) 


For arbitrary variations of the virtual solution vectors 
Sq ip) and Sk {p) (p = u,l), Eq. (68) can further be 
rearranged and put into a more compact form 


K„ 

(q < “ , ,q (<) ) 

k 12 i 


C WT 

0 

K 

(q w ,q (0 ) 


(q’-V") 

0 

£(«)T 


c ( o 


0 

0 

0 


0 


C (B) 

0 

0 


f^(0] 



q 


P 

„(«) 



q 

> = < 

p 


fit) 

l c 

p u \ 


f(u) 
^ v c 


where 


(70) 


K,j (q <u) , q (e) ) = K (<) (q <0 ) + K<? (q (!,) , q (<) ) (71a) 


K 22 (q (u) ,q w ) = K (u) (q (B) ) + K ( 2 “ ) (q (u) ,q (<) ) (71b) 

K 12 (q ( “ ) ,q (<) ) = -K^ ) (q ( “ ) ,q (<) ) (71c) 

K 21 (q w ,q w ) = -K«(qW,q») 

T , 7 ' (7 Id) 

= k ^ 2 (q “ >q f ) 

Note that the submatrices K JJ (q < “ ) ,q << ’ ) ) and 
K A<tAq w ) are dependent on the unknown 
variables, q (M) ,q (0 rendering the governing equations 
nonlinear, and that these submatrices are non- 
symmetric. The solution to this equation requires a 
nonlinear iterative solution technique that utilizes 
LU decomposition. Therefore, the Newton-Raphson 
iteration method in conjunction with Broy den’s 
automatic Jacobian matrix update procedure is 
employed. The iterative solution procedure begins with 
the initial guess of the incremental unknown variables. 
The initial guess is obtained from the linearized 
equilibrium equations, and the converged solution is 
achieved through incremental corrections to the initial 
guess. 


Numerical Results 

The present approach is first validated against the 
nonlinear finite element solution of isotropic single-lap 
joints with linear and bilinear adhesive behaviors. As 
described in Fig. 1, the adherends have identical planar 
geometries, with the width and length dimensions 
specified as W (u) = W (l) = W = 20 mm and L {u) = L {1) = 
L - 60 mm, respectively. The overlap length of the 
joint is given as L (a) = 20 mm. The thicknesses of the 
adherends and the adhesive are specified as 
h (u) = /z (0 = h = 1.5 mm and h (a) - 0.2 mm, 
respectively. These dimensions are the same as those 
considered previously by Edlund and Klarbring 9 . 

The adherends are made of aluminum with Young’s 
modulus and Poisson’s ratio specified as E {u) = E in = 
E = 70000 MPa and v (u) = v {£) = v = 0.3 , respectively. 
Also, the adhesive exhibits either a linear or a bilinear 
elastic material behavior. In the case of linearly elastic 
material behavior, the shear modulus and Poisson’s 
ratio of the adhesive are specified as G (a) = 3, 000 MPa 

and v (a) = 0.3 , respectively. In the case of bilinear 
material behavior, the material parameters of the 
adhesive are defined as G[ a) = 3, 000 MPa , G { 2 a) = 
1,500 MPa , and v (a) = 0.3 . The single-lap joint 
simply supported along the loaded ends is subjected to 
uniform tension of N 0 = 150 N/mm in 10 equal load 
increments. 
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The present analysis results are compared against 
predictions obtained from a two-dimensional nonlinear 
finite element analysis (FEA) using ANSYS, a 
commercially available program. The scaled 
deformation of the bonded single-lap joint with linearly 
elastic adhesive behavior at the final load step is shown 
in Fig. 4. 



Fig. 4 The scaled deformation of the bonded lap 
joint made of isotropic adherends with linear adhe- 
sive behavior at the last load step (N/N 0 = 1). 

As observed, the bonded single-lap joint bends 
asymmetrically as it is stretched from the right edge of 
the lower adherend. The asymmetric deformation 
occurs primarily due to the presence of eccentric 
loading, boundary conditions, and the geometrical 
coupling between the adherends. 

The transverse normal (peeling), <j zz , and transverse 
shear stresses, <j xz and o yz , in the adhesive for the case 

of linearly elastic adhesive material behavior corre- 
sponding to the final load step is illustrated in Fig. 5. 
The steep variations of the peeling , cr zz , and shearing 

stresses, cr xz and cr yz , near the edges and comers of the 

adhesive are captured successfully by the present 
approach. Note that both peeling stress cr zz and 
shearing stress a xz are symmetric along the horizontal 
and vertical centerlines of the adhesive whereas the 
shearing stress cr yz is asymmetrically distributed. 

A comparison of the peeling stress cr zz and shearing 
stress cr xz evaluated along the horizontal centerline 
(i.e., along y = 10 mm) from the present analysis with 
those of the FEA is shown in Fig. 6. The comparison 
indicates close agreement for both the linearly elastic 
and bilinear adhesive material behaviors. The present 
solution method captures the expected steep variation 
of the shearing and peeling stresses along the adhesive 
edges. 

A comparison of the slope variation at the center of 
the overlap as the load is increased for the present solu- 
tion method with that of the nonlinear FEA with 
ANSYS indicates close agreement, as shown in Fig. 7. 



Peeling stress, s 



Shearing stress, s 



Fig. 5 Variation of the stress components in the 
adhesive: peeling stress, s a ; shearing stress, s xz ; 

shearing stress, s . 


The small difference can be attributed to the modeling 
differences and the fact that FEA includes transverse 
shearing deformations in the adherends, as well as in 
the adhesive, whereas the present approach is based on 
the Kirchhoff plate theory, which excludes the 
transverse shear deformations in the adherends. 
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Normalized 
peeling and 
shearing 
stresses 




Linear adhesive 



Distance, jc, along overlap length (^=10 mm) 


Bilinear adhesive. 


Fig. 6 Comparison of the peeling stress, s a , and 
shearing stress, S xz , along the horizontal centerline, 

between the present approach and the finite element 
analysis for: linear adhesive and bilinear adhesive. 

The capability of the present approach is 
demonstrated by considering a bonded composite 
single-lap joint of angle-ply laminates (adherends) as 
shown in Fig. 1. The adherend length and width and 
the overlap length, as well as the loading and boundary 
conditions, are identical to those of the validation case 
except for the presence of tapered adherend edges. The 
lower and upper laminate edges are tapered (beveled) 
toward the adhesive edges, as shown in Fig. 1, where 
A (m) and A (0 indicate the taper lengths of the upper 
and lower adherends, respectively. 
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Fig. 7 Comparison of the rotation at the center of 
the overlap as a function of applied load between the 
present approach and the finite element analysis. 

Both adherends are symmetrically laminated and 
their angle-ply stacking sequence is given by \0 i-0 \ s , 
where 6 is referred to as the angle-ply laminate 
parameter. Each ply made of Graphite-Epoxy has 
properties of E L = 127.56 GPa, E T =11.31 GPa , 
v LT - 0.3 , and G LT - 6 GPa, with equal nominal ply 
thicknesses of t k = 0.0762 mm. Hence, the total 
thickness of the untapered laminate becomes 
2/z ( m) -2 h {l) -2h=\.2\92 mm. The thickness of the 
adhesive is specified as 2h (a) =0.12 mm. The linearly 
elastic behavior of the adhesive is defined by a shear 
modulus of G[ a) =G { 2 a) =G {a) = 0.4147 GPa, and the 
bilinear adhesive material behavior is defined by the 
parameters G[ a) = 0.4147 GPa and G[ a) = 300 GPa 

with a characteristic shear strain of y ( c a) - 0.03 . 

The effects of the tapered adherend geometry and the 
angle-ply laminate parameter, 0 , on the geometrically 
nonlinear behavior of the bonded lap joint are 
investigated by (1) varying the taper lengths, A (M) and 
A (/) , from 0 to 5 mm in five equal increments while 
specifying the value of 6 - 45 degrees (i.e, 
[+45 /- 45] 4s laminate), and (2) by varying the angle-ply 
laminate parameter, 6, as 0, 5, 10, 15, 30 and 45 
degrees while assuming an untapered adherend edge, 
i.e.,A (M) =A (/) =0. 

The effects of taper lengths on the transverse normal 
stress, <j zz , and the transverse shear stress, g xz , both 
evaluated along the horizontal centerline y = 10 mm, 
are depicted in Figs. 8 and 9, respectively, for both 
linearly elastic and bilinear adhesive material 
behaviors. As shown in these figures, the peeling stress, 
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Linearly elastic adhesive 



Bilinear elastic adhesive 


Fig. 8 Variation of the peeling stress, s a , along the 

horizontal centerline for varying adherend taper 
lengths: linearly elastic and bilinear adhesive 
behavior. 

cr zz , reduces with increasing taper length, and the 
shearing stress, a xz , reduces slightly along the 
horizontal centerline of the adhesive. The transverse 
shear stress, cr yz , evaluated along the vertical centerline 

of x = 45 mm is depicted in Fig. 10, also for both 
linearly elastic and bilinear adhesive material 

behaviors. The shear stress component near the comers 
of the adhesive edges decreases significantly with 
increasing taper length. Although not shown here, the 
transverse shearing stress, cr xz , also reduces 
significantly at the comers of the adhesive edges as a 
result of increasing taper length. As observed in Figs. 8- 


10, the bilinear adhesive material behavior yields stress 
variations identical to those of the linear adhesive 
behavior, but with relatively lower transverse normal 
and shear stress variations than those of the linearly 
elastic adhesive material. 



Linearly elastic adhesive 



Bilinear elastic adhesive 


Fig. 9 Variation of the shearing stress, s xz , along 

the horizontal centerline for varying adherend taper 
lengths: linearly elastic and bilinear adhesive 
behavior. 

Variations of the peeling stress, cr zz , and transverse 
shear stresses, <r xz and <j yz , at the final load step, 

A 0 /A“ ax =1.0, in the adhesive with linearly elastic 
adhesive material in the absence of tapered adherend 
edges are illustrated in Fig. 11. As expected, both the 
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Linearly elastic adhesive 



Bilinear elastic adhesive 


Fig. 10 Variation of the shearing stress, s yz , along 

the vertical centerline for varying adherend taper 
lengths: linearly elastic and bilinear adhesive 
behavior. 

transverse shear and peeling stresses increase rapidly 
near the edges of the adhesive. As observed in Figure 
11, all of the adhesive stress variations are asymmetric 
with respect to the horizontal and vertical centerlines. 

The adherends are symmetric angle-ply laminates 
with no material coupling between stretching and bend- 
ing that may cause twisting in the adherends and result 
in asymmetric stress distribution in the adhesive. How- 
ever, the coupling between the bending and twisting 
deformations disturbs the symmetry in the adhesive 
stresses. 



Fig. 11 Variation of the stress components in the 
adhesive: peeling stress, s a , shearing stress, s xz 

along the horizontal vertical centerline, and 
shearing stress, s yz for untapered adherends of 

[±45] 4s angle-ply laminates with linear elastic 
adhesive behavior. 


15 

American Institute of Aeronautics and Astronautics 


In order to understand this behavior, the adhesive 
stresses are evaluated along the left edge of the adhe- 
sive (i.e., x = 45) for varying angle-ply laminate 
parameter, 0 , as shown in Fig. 12. 



Peeling stress 



Shearing stress 


Fig. 12 Variation of the stress components in the 
adhesive: peeling stress, s a and shearing stress, 

s xz along the left edge of the adhesive for varying 
values of angle-ply parameter, q . 

These results are obtained for untapered adherends and 
linearly elastic adhesive behavior. As shown in Fig. 
12a, the peeling stress, cr zz , is symmetric along the 
adhesive edge for 6 = 0. As 6 is gradually increased 
from 0 to 15 degrees, the peeling stress near the upper 
corner of the left adhesive edge (i.e., x = 45 and y = 20) 


becomes higher than that near the lower comer (i.e., x = 
45 andy = 0). As 6 is further increased, the difference 
in the peeling stress near the corners of the adhesive 
edge becomes more pronounced, as shown in Fig. 12a. 
A similar behavior is also observed in the transverse 
shearing stress, cr xz , where the magnitude of stress at 

one comer is different than that at the other comer, as 
shown in Fig. 12b. Both the peeling and shearing 
stresses increase with increasing 6 . This is primarily 
due to the reduction of adherend in-plane and bending 
stiffnesses in the longitudinal direction. The lower the 
in-plane and bending stiffnesses of the adherends the 
higher the edge effects. 

Conclusions 

A semi-analytical solution procedure was presented 
for the geometrically nonlinear three-dimensional 
analysis of a single-lap joint made of composite tapered 
adherends bonded by a linearly elastic or bilinear 
adhesive. The nonlinear equations of equilibrium were 
derived based on the principle of virtual displacements 
in conjunction with the von Karman nonlinear plate 
theory for the adherends and the shear-lag theory for 
the adhesive. The bilinear adhesive material behavior 
was incorporated by computing the material parameters 
from an effective shear stress-shear strain relationship. 
The displacement fields were assumed in the form of a 
double series containing fifth-order B-spline functions 
(displacement modes) in each direction. The resulting 
nonlinear equations of equilibrium were then solved 
numerically by employing the Newton-Raphson 
incremental iterative procedure along with Broyden’s 
automatic Jacobian matrix update. 

A comparison against a two-dimensional nonlinear 
FEA solution established the capability of the present 
approach to accurately capture the steep variations of 
both peeling and shearing stresses in the vicinity of the 
adhesive edges, as well as at the corners. Furthermore, 
the variation of overlap rotations indicates that the 
stress-induced effects due to geometric nonlinearity 
were captured by preserving the nonlinear terms in the 
strain displacement relations. 

In the case of a bonded single-lap joint of angle-ply 
laminates, the tapered edges led to a considerable re- 
duction of the peeling stresses and a slight reduction of 
the shearing stress component in the longitudinal 
direction ( <j xz ). However, the increase in taper length 
reduced all the adhesive stress components at the 
comers of the adhesive region. Furthermore, the 
increase in the angle-ply laminate parameter not only 
increased the stress concentration near the adhesive 
edges but also changed the symmetric distribution of 
peeling and longitudinal (transverse) shearing stresses 
to asymmetric due to the presence of material coupling 
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between bending and twisting deformations of the 
angle-ply laminates. 

This special-purpose three-dimensional analysis 
method accurately predicts the geometrically nonlinear 
behavior of a bonded single-lap, tapered, composite 
joint subjected to uniaxial tension while accounting for 
the presence of tapering of the adherend thickness, 
changes in the overlap length of the adherends, and the 
linear and bilinear elastic behaviors of the adhesive. 


Appendix 

The B-spline functions, f (a ; t ( p , K) , with a = x,y 
and p-i,u, employed in Eq. (8) are defined 
recursively in the form 13 


( a -t (p) ) 

T (crt ip) K) = — - — - — T (or t {p) 7f-l) 

* V“’ l « ’*) (Ap) Ap )\ 1 i \ U ^cc ^ V 


/>> ~t ip) Y 

\ L a(i+K- 1 ) L a(i) ) 

(t {p) -a) 

\ L a(i+K ) U ) 

(t iP) ~t (P) ) 

ya(i+K) l a(i+ 1 ) ) 


(72) 


T i+l (a;g\K-l) 


with / = 1, 2, . . . , M {p) and k > 1 . The variable 
represents the components of the knot vector, t ip) , for 
the (K - l) st -degree B-spline functions. The knot 
vector, t {p) , is defined in terms of the Cartesian 
coordinate ( a - x, y) of the selected points as 


pp) — (Ap) Ap) 

Ap) Ap) 

' 9 a(N { P) ’ a(N^ a d) +l) ’ 

Ap) 

a(Ni d) + 2)’’ 

t ip) 1 

a(N ( a d) +2K) } 


fy(p) ,y(P) fy(P) 
..,CvCq ,Ll 2 •> 

a ip) a {p) 
V)’V41 

a {p) a {p) 


(73) 


where N {p) +1 (a - x,y;p - u,£) denotes the number 
of distinct knot points selected along the a direction in 
the lower and upper adherends. Based on the definition 
of the knot vector , t ( p , the end points (i.e., a { a p) and 
a ( p ) are repeated K times. In Eqs. (72) and (73), the 

relationship between the number of knot points and the 
extent of the B-spline functions, M ( a P) , is given as 


N { a p) =M { a p) -K + 3 (74) 

In this study, the B-spline functions with K - 5 are 
chosen to have 20 and 7 knot points in the x- and y- 


directions, respectively. The knot vectors are defined 
as 

t ( x ° = {45,45,45,45,45,46,47,49,52,55,58, 

61.63.64.64.5.65.66.69.74.80.87.5, 
95,101,106,10,110,110,110,110} 

t (u) = {0,0,0,0,0,4,9,15,22.5,30,36,41,44, 

x 1 (75) 

45.46.47.49.52.55.58.61.63.64.64.5, 
65,65,65,65,65} 

t (0 =t (u) = {0,0,0,0,0,3.33,6.66,10,13.33, 
16.66,20,20,20,20,20} 

This results in a B-spline series representation of the 
displacement components (i.e., u [ p) , u {p) , and u [ p) , 

with p = /, u) with the series terminating at M[ p) - 23 
and M {p) - 9 . 
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